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1. Introduction 

The remarkable achievement of ultraslow light^^ in a Bose-Einstein condensate (BEC), 
using electromagnetically induced transparency^^ (EIT), inspired many intriguing and at- 
tractive applications^"^, in particular quantum optical dynamic memories^i^ 1 ^^ 1 ^. A 
key component of a quantum information processor would be the dynamic memory storage 
device. BECs have the potential to be used in such an application. Consider a sequence of 
information bits carried by a train of optical pulses. Upon injection into a BEC, the coherent 
optical information can be stored for a period given by the transit time of the pulses. In the 
BEC phase, the demonstrated ultra-slow group velocities allow for long storage periods of 
the order of few microseconds. Operational performance of such a device can be character- 
ized by the overall bit storage capacity of the BEC, which can be defined as the number of 
optical pulses that can be simultaneously stored during the storage time in the BEC. 

In order to make quantum memories more practical for quantum information processing, 
it is useful to enhance their bit storage capacity. This would be possible, in principle, by 
using pulses temporally as short as possible so that more pulses can be injected into the 
BEC during the storage time^. Early ultraslow light experiments use pulses of widths few 
microseconds. Shorter pulses would have longer frequency bandwidths, beyond the EIT 
window. EIT window can be modified by an external magnetic fields. Very recently, it has 
been experimentally demonstrated that EIT in a resonant atomic medium is possible for a 
comb of picosecond pulses^. Promising schemes to introduce large controllable time delays 
for such ultrashort pulses has been proposed where EIT window is overcome by using spatial 
and temporal processing 19 . Distortions in the shape of a large-bandwidth pulse have been 
experimentally studied in warm Rb vapor—. Shapes of the stored pulses will depend on 
the dispersive characteristics of the BEC. In particular, each pulse in the train will undergo 
dispersive broadening. The temporal separation of the pulses should be adjusted in such 
a way that after broadening in the BEC, consecutive pulses do not have significant spatial 
overlap. This sets the limit to the bit storage capacity of the BEC and can be estimated from 
the overall group delay and pulse broadening. It is therefore necessary to investigate the 
dispersive effects inside a BEC and develop accurate models that describe pulse propagation. 

In this paper, we consider propagation of short pulses with widths in the microseconds 
to nanoseconds range. In particular, we investigate the role of dispersion on the temporal 
characteristics of optical pulses propagating inside a BEC. Analytical formulas are first 
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summarized for the constitutive relations from which the group delay, absorption coefficient, 
and the second-order group delay dispersion can be calculated. We find that third and 
higher-order dispersion effects are negligible. The influence of the spatial inhomogeneity of 
the BEC, the interatomic scattering interactions and the temperature on the pulse shape 
are taken into account in our treatment. We show that just below the critical temperature 
of condensation, broadening becomes maximum. Finally, we give analytical formulas for the 
bit storage capacity of a BEC and determine its optimum performance. 

2. Propagation and dispersion of slow-light pulses in an atomic ensemble under 
EIT scheme 

We consider a gas of N three-level atoms interacting with two laser beams in A configuration 
as shown in Fig{TJ In this scheme, it is assumed that the only dipole forbidden transition is 
the one between the lower levels. The upper level is coupled to the lower levels via a strong 
drive field with frequency uj c and a weak probe field of frequency u p . This particular ar- 
rangement leads to the well-known effect of electromagnetically induced transparency (EIT) 
in which the absorption of the probe field at resonance is cancelled. The underlying physical 
mechanism of EIT is the quantum interference induced by the drive field. Quantum coher- 
ence established between the lower level doublet put them in a dark state where excitations 
to the upper level become impossible. The transition probabilities destructively interfere 
for the two possible excitation routes to the upper level. In addition to turning an optically 
thick, opaque, medium transparent, EIT has recently been used for achieving ultra-slow light 
velocities, owing to the steep dispersion of the EIT susceptibility near the probe resonance. 
Susceptibility x f° r the probe transition can be calculated as a linear response as most of the 
atoms remain in the lowest state. Assuming local density approximation, neglecting local 
field, multiple scattering and quantum corrections and employing steady state analysis, \ is 
found to be£ 

pKI 2 ^A + r 2 /2) 

X " e h [(T 2 /2 + iA)(r 3 /2 + iA) + n2/4]' u 
where A = cg> 31 — uj p is the frequency detuning of the probe field with frequency u p from the 
resonant electronic transition u) 3 i- For the cold gases considered in this paper and assuming 
co-propagating laser beams, Doppler shift in the detuning is neglected. It is assumed that 
the drive field is at resonance. Q c is the Rabi frequency of the drive field; //31 is the dipole 
matrix element between states |3) and |1) which can also be expressed in terms of resonant 
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wavelength A31 of the probe transition via fi\ x = SeohX^^ /8ir 2 with 7 is the radiation decay 
rate between |3) and T2 and T3 denote the dephasing rates of the atomic coherences 
between the appropriate states, p stands for a given atomic density. 

Probe beam propagates along the condensate axis in the z direction. We treat the axial 
propagation of probe beam under paraxial approximation where the paraxial effects and 
diffraction losses are ignored. This approximation is valid when the probe beam radius is 
much larger than the radial size of the atomic cloud. In practice, using a pinhole and a 
flipper mirror, only that portion of the probe beam passing through the thin central column 
region of the cloud is selectively monitorec—. To develop a one dimensional wave equation 
including the high order diffraction effects in the dispersive medium for the short probe 
pulse, we follow Ref.— . 

The probe pulses we shall consider in this paper have carrier frequencies in the order of 
10 15 Hz, while their spectral widths are less than of the order of 10 9 Hz. We will study their 
propagation under slowly varying phase and envelope approximations. The wave equation 
that would govern the propagation of a short pulse in a dispersive medium can be found to 
be^i 

BE 1 dE d 2 E 

where a determines the pulse attenuation; v g is the group velocity, and b 2 is the second-order 
dispersion coefficient, or the group velocity dispersion. We have also calculated the third 
order dispersive coefficient 63 and found it to be 7 orders of magnitude less than 62 for typical 
experimental parameters^. The significant coefficients for the short pulse propagation can 
be calculated from the susceptibility using the relations^ 

%ix 

« = -txW, (3) 
v g c A Duj K ) 
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Here we take A31 = A, u p = u and u 3 i = ujq for convenient notation. It may be noted 
that these expressions lead to complex parameters in general, as the x is a complex valued 
function and the parameters are merely mathematical coefficients in Eqj2j At EIT resonance 
however they lead to physically meaningful, well-defined, real valued absorption coefficient, 
group velocity and dispersion coefficient. EIT susceptibility exhibits steep normal dispersive 



behavior around the resonance, which allows for the second term in EqJUbecome much larger 
than the first term. This leads to substantial reduction in the group velocity of the pulse. 
Such a slow pulse propagates through the medium without much absorption due to the small 
imaginary EIT susceptibility at resonance. On the other hand, for short pulses the pulse 
shape would also be influenced by the higher order dispersive coefficients, starting from the 
group velocity dispersion. Different frequency components would have different velocities 
and broadening of the pulse may become significant. Let us emphasize that essential physical 
foundations of these effects are not related to the density p but the dispersive properties 
of the EIT susceptibility. Indeed, slow laser pulses have been observed in various media, 
including hot rubidium vapor— >22, ultracold Bose gas^, and in solid crystals 2 ^. When Rabi 
frequency for the coupling field is sufficiently large such that fl c 3> r 2) 3, we can calculate 
these parameters characterizing propagation and dispersion of slow EIT pulses explicitly as 
follows 



a = e nxni ' (6) 

ce hQ 2 

^9 r> i 19 ' 

2a>3i|/i 31 | 2 p 
.8nT 3 \p 31 \ 2 p 
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It should be noted that the result for v g is valid when 



r 2 , 3 < < y/A7rc\n3i\ 2 p/e h = y ^cA 3 7 p. (9) 

For p ~ 10 20 — 10 21 m~ 3 , the upper limit would be 57 — 157. In this paper we consider Q c = 
O.567 — 57. Corresponding electric field of the coupling laser then would be in the range 10 2 — 
10 3 V/m, with power densities 13 — 1300 W/m 2 . Q c = O.567 is used in typical experiments^. 
We shall discuss possible enhancement of bit storage capacity by selecting particular fl c 
within this range. According to their dependencies on the Q c , we see that increasing fl c 
reduces the dispersive effects significantly, but simultaneously it reduces the storage time of 
the pulses in the condensate as the group velocity would rise. It is illuminating to realize 
that decrease of delay time (storage time) is much slower than the decrease of dispersion 
coefficient with the Q c . This can be exploited to find a critical fl c for which number of probe 
pulses injected into the medium within the storage time would be optimized. 
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We note that for extremely large Q c , one should use the following expression instead 



For a uniform density medium the Eqj2]can be solved analytically 2 ^. An initial temporally 
short gaussian pulse after propagating in the medium of length L is then found to be delayed 
with respect to a reference pulse propagating in vacuum by t d = L/v g — L/c. During that 
time, the pulse would broaden due to second-order dispersion. Leaving the medium, it would 
grow to a final width of 



with t is the initial temporal width of the pulse, and zq = — vttq/62- For L ^> zo we get 
t(L) = \b 2 \L/nTo. 

In general, p = p(r, z) is a spatially varying function, whose profile depends on tempera- 
ture, inter-atomic interactions and confinement potential. Assuming slow spatial variations 
within an optical wavelength the wave equation remains the same. The only effect of spa- 
tial inhomogeneity would be to give group velocity and its dispersion local character. As a 
result, the time delay td can be calculated through a spatial averaging of the group velocity 
field^. Experimentally measured group velocity is then operationally defined by v g = L/td, 
where the effective axial length of the medium is evaluated by 



The axial length L is an effective length corresponding to the axial width of the density 
distribution. It may be noted that more exact treatments can be used to determine the 
effective length of an interacting BEC as discussed in detail very recently 27 . There is no 
simple operational definition for the broadening of the pulse. One has to either determine the 
output pulse shape from exact numerical simulations of the wave equation or refer to uniform 
density results for qualitative studies. While it is the atomic density which determines 
the local group velocity and its dispersion, after spatial averaging, the group velocity and 
broadening of the pulse would become explicitly dependent on the temperature, interactions 
between the atoms as well as the other physical parameters due to the confinement potential. 
In order to elucidate these dependencies, it is necessary to define p concretely and discuss 
how it is shaped by such physical parameters. We shall specifically consider an atomic 
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Bose-Einstein condensate as the medium of propagation in this paper and review its density 
profile in the next section. 

3. Density profile of a trapped semi-ideal Bose gas at finite temperature 

At low temperatures a Bose gas can be considered to be composed of two components. The 
first component becomes the condensate part and the other is a thermal gas background. In a 
harmonic trap, spatial overlap between the components becomes small, which makes effects 
of collective elementary excitations weaker relative to those in homogeneous gases. This 
allows for semiclassical descriptions of density distributions in terms of the self-consistent 
Hartree-Fock model, which simplifies more general Hartree-Fock-Popov description^. When 
the thermal component is dilute enough, one may further neglect the atomic interactions in 
the thermal component. Taking into account mean field repulsion only among the condensed 
atoms, an analytical explicit description of a partly condensed gas was developed and denoted 
as the semi-ideal model^S. We will now briefly review this model, where condensate density 
is evaluated under Thomas-Fermi approximation^^ to Gross-Pitaevskii equation^ 2 -^, and 
the thermal gas density is calculated semi-classically 2 ^ 1 ^. The total density at a temperature 
T is then written to be 



chemical potential; 6(.) is the Heaviside step function; g n (x) = £j x J /j n is the Bose function; 
At is the thermal de Broglie wavelength; (3 = l/k B T; z = exp is the fugacity, and Tq 
is the critical temperature. The external trapping potential is V(r) = (m/2)(uj 2 r 2 + u 2 z 2 ) 
with to r the radial trap frequency and u> z the angular frequency in the z direction. \i is 
determined from N — j d 3 fp(r). At temperatures below T c this yields 29 



where htf is the chemical potential evaluated under Thomas- Fermi approximation and the 
condensate fraction is given by 




(13) 



where Uq = 4irh 2 a s /m; m is atomic mass; a s is the atomic s-wave scattering length; ji is the 
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with x = T/T c , and ( is the Riemann-Zeta function. The scaling parameter s, characterizing 
the strength of atomic interactions within the condensate, is calculated to 

be 28^9 



Here, ah = \jh/m{uo z uj'^) 1 / 2 ' denotes the average harmonic oscillator length scale. At tem- 
peratures above the Tc, p can be determined from Li^(z) = ((3)/x 3 , where Li^(.) is the 
third-order polylogarithm function. The semi-ideal model has a wide-range of validity in 
representing density distribution of a trapped Bose gas at finite temperature provided that 
s < 0.4M. At the same time the interactions are assumed to be strong enough to en- 
sure p ^> hu r)Z so that kinetic energy of the condensate can be neglected according to the 
Thomas-Fermi approximation. In typical slow-light experiments in cold atomic gases, s re- 
mains within these limits. We shall use the semi-ideal model given density profile p in our 
investigations of propagation and dispersion of slow pulses and present our explicit results 
in the next section. 

4. Results and discussions 

In our numerical calculations of the theoretical results presented in the preceding section, 
we shall specifically consider a gas of N = 8.3 x 10 6 23 Na atoms for which M = 23 amu, 
A = 589 nm, 7 = 2vr x 10.01 MHz, T 3 = O.57, T 2 = 2?r x 10 3 Hz, and a s = 2.75 nm. For the 
parameters of the trapping potential, we take u r = 2tt x 69 Hz and u z = 2tt x 21 Hz as in 
Ref. 1 . The coupling field Rabi frequency is taken to be Q c = O.567 1 . Critical temperature 
for Bose-Einstein condensation of such a gas is found to be Tc = 424 nK. We illustrate the 
spatial density profile of the condensate in the axial (z) direction at T = 43 nK in Figj2j 

Due to spatially inhomogeneous density, the linear dielectric susceptibility would also 
be spatially inhomogeneous. In this case, group velocity defined by Eq.(J7j) would have a 
local character. Under EIT conditions, when the light pulse enters the condensate from the 
thermal component of the gas, its group speed exhibits a dramatic slowing down as shown 
in Figj3] Here we consider resonant probe pulse with A = 0. Within the condensate region, 
at such low temperatures, the group velocity remains roughly at the same ultraslow value. 
The pulse rapidly accelerates to high speeds when it leaves the condensate at the interface 
to thermal component. Before and after the thermal component, the pulse is assumed to 
be propagating in vacuum. In practice, group velocity is measured in terms of time delay 
of the pulse with respect to a reference pulse propagating in vacuum over the same distance 
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with the gas. To make comparisons with this operational definition of the group velocity, 
one needs to make careful spatial averaging of the theoretical group velocity 26 . 

The other coefficients in the wave equation are also locally defined. The coefficient of 
absorption and the second-order dispersion coefficient are shown in FigfJ] and in FigfJ3, 
respectively. Under EIT conditions, imaginary part x" °f X is small, despite the large 
condensate density. For large Q c we have x" ~ P^2/^ 2 - For short pulse propagation we 
shall see that the major source of reduction in the peak of the pulse would be the temporal 
broadening due to dispersion. The absorption over the small size of the BEC is negligibly 
small. Taking a = 10 3 m _1 and L = 100 jttm, we can roughly estimate attenuation of the 
pulse by exp (— otL), which yields about 90% transmission. The loss term can be made 
smaller by using larger Rabi frequencies for the coupling control field. The positive sign of 
a reflects its effect on the pulse as the decrease of the pulse intensity. The broadening of 
the pulse on the other hand is independent of the sign of b 2 , as can be deduced from EqJTTl 

The region over which dispersive effects are most strong is a small region near the center 
of the cloud as can be seen in FigfSl The size of this region can be determined by the width 
L of the density distribution, since 62 is proportional to p. It is quite smaller than the 
classical boundary ^J2p/mu 2 . As an example, L « 26 /iin in FigfSl while the condensate 
is extended over a size of ~ 140 pm. In such a small region, spatial variations of p can be 
ignored so that p « p(r — 0, z — 0) = po- When a pulse of width tq enters into this strongly 
dispersive zone, its final width would approximately broaden to t(L) = T ^Jl + n 2 (p L) 2 /tq, 
where k = 3r 3 7A 2 /7r 2 f^. Broadening of the pulse depends on the density profile of the 
medium through its characteristic parameters, its peak value po, and its width L. For small 
dispersion, change in the pulse width would be proportional to (pqL) 2 . For the given p within 
the semi-ideal model, L demonstrates the well-known temperature variation 2 ^ as shown in 
Figj6j At a given temperature, L has been evaluated by Eqfl2l At very high temperatures, 
T ^> Tc, classical thermal radius would be ~ \/T according to the equipartition theorem. 
In the temperature range of FigJHl beyond Tc, L is changing linearly with T. It might be 
illuminating to rewrite EqJT2j as L = J (4tt/N)((z 2 )c + (z 2 )t)- Here, (z 2 )c,t stands for a 
mean square axial distance evaluated using either condensate (C) or thermal density (T) 
distribution. The sharp drop of the width of the distribution just below T c is due to the 
emergence of condensed component with relatively high density about the center of the 
cloud. As T becomes lower, thermal component shrinks, while the condensate part expands 
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slowly about the center of the cloud. L eventually saturates down to Thomas-Fermi width at 
zero temperature. The thermal behavior of the peak density is also well-known2£>22. At high 
temperatures it shows ~ 1/T 3 / 2 decay consistent with classical gas behavior. Just below To, 
Po grows rapidly due to emerging relatively dense condensate about the center of the cloud. 
At extremely low temperatures, po saturates to Thomas-Fermi density. 

Figd describes dependence of broadening on temperature and interatomic interaction. 
At extremely low temperatures, L and po are determined by the condensate part and they 
change slowly with T. For increasing T, thermal part expands beyond the boundaries of 
the condensate, so that L will be contributed significantly by the thermal part. Except 
at temperatures close to To, Po is dominated by the condensate part, which exhibits slow 
dependence on T in comparison to that of L. As a result, the product p L grows, following 
the expansion of thermal component, up to a point just below Tc, where po dramatically 
drops to dilute thermal density value due to rapid disappearance of the dense condensate 
component. After that peak point, po shows faster variation with T relative to that of 
L. Hence, the shape of broadening curve is characterized by the thermal behavior of po- 
At temperatures higher than Tc, broadening decreases with temperature. At such high 
temperatures, L grows slowly (almost linear) while po decreases like 1/T 3 / 2 . Thus their 
product becomes a decreasing function with T. Summarizing, the broadening is determined 
by the density profile through its characteristic product PqL. As po and L exhibit competing 
thermal behaviors, a peak arises just before the Tc, a signature of emerging condensate being 
a strongly dispersive region. Before the peak point, it is L which shapes the temperature 
dependence of the broadening. After the peak point, broadening behaves analogous to the 
thermal behavior of po- 

Earlier studies indicate that group velocity depends on the strength of atomic interactions 
within the condensate which is characterized by the s-wave scattering length a s . As the a s 
grow longer, the group velocity of the pulse within the BEC becomes faster—. In order to 
discuss the influence of atomic interactions on broadening of the pulse, we first note that 
below Tc the atomic energies are so low that, there is only a s , s-wave scattering length, 
to characterize the strength of interactions within the condensate component. At higher 
temperatures beyond Tc, the interactions has no effect within the semi-ideal model where 
thermal gas is dilute and non-interacting. In FigJTl we see that as the a s increases, the 
broadening decreases when T < Tc- According to semi-ideal model density profile, we get 
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Pq = fi/U at T < T c , from which we conclude p ~ l/ a ^ 5 - This is intuitively expected, 
as the increase of repulsive atomic interactions results in lower densities. The width of the 
atomic cloud would slightly increase with a s . The increase of the size of the condensate 
can be estimated from the Thomas- Fermi axial radius given by ~ a l^ 5 - Thus, 

we see that the increase of L would be slower than the decrease of the p , so that p L and 
hence the broadening would reduce with a s . This conclusion remains the same at higher 
temperatures where the thermal length scale becomes dominant in determining L, and it 
is not influenced by a s . Thus, broadening decreases due to decrease of p . Beyond T c , the 
density profile is semi-classically determined and is independent of a s . The semi-ideal model 
density profile becomes not as good description of the condensate when the interactions 
characterized by the scaling parameter s becomes larger than 0.4, corresponding to a s ~ 7 
nm. For 23 Na cloud, a s = 2.75 nm. The higher a s may be achieved by utilizing Feshbach 
resonance technique 3 ^" 3 ^ 3 ^ or by considering alternative atoms. 

Our approach in determining the pulse broadening using the effective medium width 
L , as we will show below, gives qualitatively accurate behavior of the broadening but 
overestimates it. In an actual pulse propagation, a pulse of width r would be prepared far 
from the gaseous medium. Before reaching to the most dispersive zone about the center of 
the cloud, the pulse would pass through the thermal cloud and then through the edges of 
the condensate. Due to the inhomogeneous density profile, when the pulse arrives to the 
center of the cloud, it would already be much broader than r . As a result it would suffer less 
from dispersion than the one predicted above. Despite the lack of actual input value To for 
our analysis above, the overall qualitative behavior is expected to be accurate. To test this 
expectation and to justify our effective central dispersive zone treatment agree qualitatively 
with the exact pulse behavior, we now investigate the variation of the pulse width with 
distance of propagation, taking into account the full spatially inhomogeneous density profile 
of the gaseous medium. We assume a Gaussian pulse of the form exp {—(3{t — to) 2 ) & t initial 
time t , , where yffi is the pulse width, and propagate it numerically. A dimensionless 
form of Eq.flS]) is solved via finite difference Crank-Nicholson space marching scheme. The 
Crank-Nicholson scheme is less stable but more accurate than the fully implicit method; it 
takes the average between the implicit and the explicit schemes^ 9 -. We use forward difference 
scheme for the position and central difference scheme for the time. Discrete equations in 
matrix form are solved using Thomas algorithm 39 which is analogous to a fast Gaussian 
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elimination method for tridiagonal matrices. 

We have performed extensive numerical simulations at various T and confirmed that de- 
spite being an overestimation, the qualitative behavior predicted in FigfT] correctly describes 
the results obtained by the exact solution. We present typical results of our simulations in 
FigslSlTTTl For a nanosecond pulse, at T = 43 nK, broadening is ~ 3 /is which can be seen in 
Figj9j This is much smaller than the result found by the uniform density calculation, which 
gives ~ 100 fis at the same temperature in FigJTl More moderate broadening is obtained for 
a pulse of width 0.1 /is and shown in FigJTQl As a rule of thumb, we assume that broadening 
is not significantly large if it is less than a factor of 2. A microsecond pulse broadens by 
a factor of ~ 1.7 as can be seen in Fig{TTJ For longer pulses therefore we conclude that 
broadening is not significant. 

In early ultraslow light experiments microsecond pulses were used^. Previous theoretical 
models ignore dispersion effects in the pulse propagation. While they can be used to describe 
observed group velocities, they cannot explain the transmitted pulse shape adequately. High- 
order dispersion, absorption, transverse diffraction, and off-resonant transitions to upper 
levels not included in the three-level model cause imperfect transition, and result in reduction 
of the transmitted peak intensity. We have already argued that absorption, though still 
present even under EIT conditions, makes a little effect over the small BEC size. Transverse 
diffraction due to paraxial pulse propagation is also not too strong^. To characterize the 
contribution of second-order dispersion on the reduction of the peak intensity, propagation 
of a microsecond pulse is simulated and presented in Fig{TTJ This result indicates that in 
addition to off-resonant transitions to other levels, broadening of the pulse also contributes 
to the observed drop of the transmitted peak intensity. 

While microsecond pulses are good to minimize dispersive effects, shorter pulses may 
be more attractive for optical information processing. The question we would like to now 
address is if such shorter pulses can still be stored as efficiently via EIT scheme in atomic 
ultracold gases. We may characterize coherent optical information storage capacity through 
a parameter given by C = L/2v 9 t(L). It measures number of probe pulses simultaneously 
present in BEC without significant overlap. The storage time is taken to be the delay time 
of the pulses, t s = L/v g . We choose an ideal repetition rate of the pulse train to be l/2r(L). 
The pulses are assumed to be stored in the central region of the cloud, whose length L is 
given by EqJT2l It is the same with the effective central dispersive zone considered in the 
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analytical broadening calculation. In this region p fa po is approximately uniform at its peak 
value, leading to spatially homogeneous v g , that can be seen in Figj3l More explicitly, C 
becomes 

C = L (17) 

2r - v /47r 2 ^/9A4 7 2p2 + 4^2/^2,^4 

The second term under the square root in the denominator represents the second-order 
dispersive contribution. When dispersion is not too strong, C could be made higher by 
increasing p. On the other hand, it is not very simple task to increase p for BEC. Besides, it 
implicitly and indirectly affects the dispersive term via L. For a given r , more explicit and 
direct control parameter would be Q c . C exhibits a non-trivial dependence on Q c , which is 
illustrated in FigJT2"l We observe that for given atomic cloud, and r , there is an optimum 
choice of fl c . Fig{T2]also compares C for different r . We deduce that, for a wide range of 
Q c it is indeed possible to inject more pulses by using shorter r . However, reducing r from 
1 pis to 0.01 pis, one cannot hope to increase the capacity 100 times more. Instead, due to 
dispersive effects, only by using a critical coupling laser Rabi frequency, storage capacity 
can be made maximum. The maximum achievable capacity C max is independent of To, and 
is found to be 



Cmax ~i~zzrf' (18) 

which is evaluated at the critical coupling Rabi frequency, Q c q, given by 

These results suggest that shorter pulses can be stored in the cloud as efficiently as microsec- 
ond pulses. The storage time t s0 when Q c = Q c0 becomes 

(2 °> 

We note that analytical result of C underestimates the actual value of it, as t(L) is 
overestimated. Due to spatial inhomogeneity, r , pulse width just before the central zone, 
would be larger than the width of the original pulse prepared away from the atomic cloud, 
so that t(L) would be smaller. v g would remain the same under any pulse width change. As 
a result, actual C would be higher than the one predicted above. The analytical expressions 
are mainly good for predicting qualitative behavior of C with VIq, and providing a lower 
limit for its value. 
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5. Conclusion 

We have examined dispersive effects on short pulse propagation through a semi-ideal BEC 
under EIT scheme. We have found that third and higher order dispersion coefficients are 
negligibly small for current experimental systems. Distortion of the shape of a short pulse 
in the atomic cloud would be a temporal broadening arising from second-order dispersion 
per se. For pulses of widths greater than microseconds we have found that second order 
dispersion is also negligible. For microsecond pulses the broadening is about 1.7 of the 
initial width. As a rule of thumb, when the pulse length is broadened by less than a factor 
of 2, the pulse can be considered to be not significantly distorted. 

Second order dispersion of the atomic cloud depends on its density profile, which is 
determined mainly by temperature and atomic interactions. For nanosecond pulses, we have 
presented the dependence of the broadening on the temperature and interatomic scattering. 
The broadening is found to be making a peak just below the critical temperature T C - As 
the scattering length increases, the broadening decreases at all temperatures, resulting in a 
lower peak. 

Furthermore, we have numerically examined the effect of spatial variations of the atomic 
cloud. Our results indicate that due to spatial inhomogeneity the broadening is less than that 
predicted by the calculations based upon effective central dispersive zone. As the pulse is 
already broadened relative to its original width, when it arrives at the central zone, it suffers 
less from the high dispersive properties of this region. Qualitatively, we have found effective 
central zone estimations describe the behavior of the broadening very well. Our numerical 
simulations showed that the transmitted pulse shape in the ultraslow light experiments starts 
to get distorted for pulses of widths shorter than or equal to microsecond. 

Finally we have shown that despite the dispersive effects it is still possible to achieve 
significant increase of bit storage capacity by carefully choosing a critical coupling field 
Rabi frequency. Though shorter pulses allow for higher capacity for a wide range of Rabi 
frequencies of the coupling field, the maximum achievable capacity is independent of the 
initial pulse width of the probe pulse and determined by the properties of the atomic cloud 
only. 

Our results may be used for measuring T c and a s of the BEC medium by comparing the 
group speed and the broadening of the probe pulse with a reference pulse propagating in 
the absence of medium. One may also exploit Feshbach resonances to tune a s to control 
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broadening. Further uses of our results may be found in pulse shape engineering, frequency 
filtering, and enhancing capacity of dynamic optical memories. 
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List of Figure Captions 

Fig. 1. Schematic energy level diagram of a three- level atom interacting with two laser 
beams in A configuration. 

Fig. 2. Axial density profile of 23 Na the Bose-Einstein condensate of N = 8.3 x 10 6 atoms at 
T = 43 nK. p is scaled by the peak density. Trap parameters are chosen to be u r = 2n x 69 
Hz and oo z = 2ty x 21 Hz. 

Fig. 3. Spatial dependence of the local group velocity of a resonant probe pulse along the z- 
axis, with A = 0, propagating through a 23 Na the Bose-Einstein condensate of A^ = 8.3 x 10 6 
atoms at T = 43 nK under EIT scheme. The parameters used are M = 23 amu, a s = 2.75 
nm, A = 589 nm, 7 = 2tt x 10.01 MHz, T 3 = O.57, Q c = O.567, T 2 = 2vr x 10 3 Hz. 
Fig. 4. Position dependence of the absorption coefficient a along the z-axis. The parameters 
are the same with those of Figj3] 

Fig. 5. Axial spatial profile of the second-order dispersion coefficient 62 for the parameters 
same with those of Figj3l Vertical axis is scaled by the peak value of b 2 which is ~ 0.108 
s 2 /nm. 

Fig. 6. Thermal behavior of the peak p and the effective axial length L of the atomic 
density distribution, given by the semi-ideal model. The parameters are the same with 
those of Figj3l 

Fig. 7. Thermal behavior of the broadening of a nanosecond pulse obtained analytically 
by an effective dispersive zone treatment. The solid, dashed, dotted and dot-dashed lines 
are for a s = 2.75, 3.75, 5.75 nm, and a s = 7 nm, respectively. The other parameters are the 
same with those of FigJHl 

Fig. 8. Thermal behavior of the broadening of a nanosecond pulse determined by numerical 
simulations. The solid and dashed lines are for a s = 2.75 nm and a s = 7 nm, respectively. 
The other parameters are the same with those of FigJSJ 

Fig. 9. Propagation of a nanosecond pulse through the interacting BEC. Time (t) is scaled 
by 0.4 ps. and position (z) is by 10 pm. The parameters are the same with those of FigJHl 
Fig. 10. Propagation of a pulse of width 0.1 ps through the interacting BEC. Time (t) is 
scaled by 0.4 ps. and position (z) is by 10 pm. The parameters are the same with those of 
FigEl 

Fig. 11. Propagation of a microsecond pulse through the interacting BEC. Time (t) is scaled 
by 0.4 ps. and position (z) is by 10 pm. The parameters are the same with those of Figj3l 
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Fig. 12. Coherent optical information storage capacity C (dimensionless) as a function of 
coupling field Rabi frequency Q c , scaled by 7. Parameters used are L = 100 /xm, r = 1 /xs 
(solid line), r = 0.1 /xs (dashed line), r = 0.01 /xs (dot-dashed line), p = 10 21 m~ 3 . The 
other parameters are the same with those of FigJSl 
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Fig. 1. Schematic energy level diagram of a three-level atom interacting with two laser 
beams in A configuration. tarhanFl.eps. 
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Fig. 2. Axial density profile of 23 Na the Bose-Einstein condensate of N = 8.3 x 10 6 atoms at 
T = 43 nK. p is scaled by the peak density. Trap parameters are chosen to be u r = 2ir x 69 
Hz and to z = 2ir x 21 Hz. tarhanF2.eps. 
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Fig. 3. Spatial dependence of the local group velocity of a resonant probe pulse along the z- 
axis, with A = 0, propagating through a 23 Na the Bose-Einstein condensate of N — 8.3 x 10 6 
atoms at T = 43 nK under EIT scheme. The parameters used are M = 23 amu, a s = 2.75 
nm, A = 589 nm, 7 = 2tt x 10.01 MHz, T 3 = O.57, Q c = O.567, T 2 = 2tt x 10 3 Hz. 
tarhanF3.eps. 
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Fig. 4. Position dependence of the absorption coefficient a along the z-axis. The parameters 
are the same with those of FigOJ tarhanF4.eps. 
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Fig. 5. Axial spatial profile of the second-order dispersion coefficient 62 for the parameters 
same with those of Figj3l Vertical axis is scaled by the peak value of 62 which is ~ 0.108 
s 2 /nm. tarhanF5.eps. 
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Fig. 6. Thermal behavior of the peak p an d the effective axial length L of the atomic 
density distribution, given by the semi-ideal model. The parameters are the same with 
those of FigOl tarhanF6.eps. 
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Fig. 7. Thermal behavior of the broadening of a nanosecond pulse obtained analytically by 
an effective dispersive zone treatment. The solid, dashed, dotted and dot-dashed lines are 
for a s = 2.75, 3.75, 5.75 nm, and a s = 7 nm, respectively. The other parameters are the 
same with those of FigfJH tarhanF7.eps. 
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Fig. 8. Thermal behavior of the broadening of a nanosecond pulse determined by numerical 
simulations. The solid and dashed lines are for a s = 2.75 nm and a s = 7 nm, respectively. 
The other parameters are the same with those of FigJSJ tarhanF8.eps. 
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Fig. 9. Propagation of a nanosecond pulse through the interacting BEC. Time (t) is scaled 
by 0.4 /is. and position (z) is by 10 /im. The parameters are the same with those of Figj3l 
tarhanF9.eps. 
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Fig. 10. Propagation of a pulse of width 0.1 /is through the interacting BEC. Time (t) is 
scaled by 0.4 /xs. and position (z) is by 10 fim. The parameters are the same with those of 
Figj3j tarhanF10.eps. 
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Fig. 11. Propagation of a microsecond pulse through the interacting BEC. Time (t) is scaled 
by 0.4 /is. and position (z) is by 10 /im. The parameters are the same with those of Figj3l 
tarhanFll.eps. 
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Fig. 12. Coherent optical information storage capacity C (dimensionless) as a function of 
coupling field Rabi frequency Q c , scaled by 7. Parameters used are L = 100 jttm, r = 1 /is 
(solid line), r = 0.1 /xs (dashed line), r = 0.01 /is (dot-dashed line), p = 10 21 m~ 3 . The 
other parameters are the same with those of Figj3l tarhanF12.eps. 
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